The Equation of State 





The equation of state is the function that relates the pressure to 
the density, molecular weight, and temperature at any place in 
the star. Since it is a solely an internal property of the gas, it can, 
in principle, be computed once externally, and used via a lookup 
table, i.e., P = P(p, u, T). 





The net pressure can be divided into three components, pressure 
from radiation, pressure from ions, and pressure from electrons. 
The latter may not obey the ideal gas law due to the effects of 
degeneracy. 


Radiation Pressure 





Normally, the pressure due to radiation in stars is small. How- 
ever, in the centers of massive stars, the energy generation is large 
enough to make the term substantial. ‘To compute the radia- 
tion pressure, recall that from kinetic theory, the pressure of an 
isotropic gas comes from the transfer of momentum, i.e., 
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Since the energy density of a photon gas is 
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then the number density of photons, n = u,/hv, in terms of the 
momentum, p = hv /c, is 
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The pressure from the photon gas is therefore 
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or, letting x = pc/kT (i.e., p= kT x/c), 
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The integral in (7.1.2) can be evaluated analytically (although I 
wouldn’t recommend it), and is 
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Thus 


where a is the radiation constant (a = 40/c = 7.56471 x 107*° 
ergs cm~? °K~*), 


Ion Pressure 





Under almost all circumstances, the pressure from ions can be 
evaluated through the ideal gas law 
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The exceptions are 


Ion Crystallization. At extremely high densities and low tempera- 
tures, Coulomb interactions between ions can become important. 
The ions are then constrained to move within a lattice, and their 
properties change. This should happen when the thermal energy 
of each ion becomes less than the inter-ion Coulomb energy, t.e., 
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If we note that the volume per ion, V = 4/3 mr3,, = 1/nion, then 
the ratio of Coulomb energy to thermal energy is 
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(Actually, crystallization begins to become important when re ~ 
100.) The melting temperature is thus 
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Except for evolved white dwarfs, stellar temperatures are always 
greater than the melting temperature. 





Ion Degeneracy. Because the mass of an ion is < 2000 times that 
of an electron, the momentum vectors associated with ions are 
much larger. (In the case of protons, by a factor of \/m,/me.) 
Thus, the phase space available for ions is much greater than for 
electrons, and ion degeneracy is not a problem. Note also that 
a-particles are not Fermions, and thus degeneracy rules do not 
apply to compact helium cores. 


Electron Pressure 





For main sequence stars, the pressure from electrons obeys the 
ideal gas law. But for more evolved stars, electron degeneracy may 
be important. To understand how and why degeneracy works, first 
compare the Maxwellian velocity distribution 
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with the Pauli exclusion principle. According to Pauli, each elec- 
tron occupies a region of phase space 
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and only two electrons (spin-up, spin-down) can occupy the same 
region of phase space at the same time. Therefore, since the vol- 
ume of phase space available to electrons with momenta between 
p and p+ dp is 4rp?dp, the maximum phase-space density of these 
electrons is 

(4rp*) 
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For high densities, the number of electrons at a given momentum 
state may exceed the exclusion condition. ‘The electrons are then 
forced into higher momentum states, until the Fermi momentum, 
pf is reached, t.e., 
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This equation defines the Fermi momentum. It is the minimum 
momentum required to fit ne dV electrons into their lowest states. 


COMPLETELY DEGENERATE, NON-RELATIVISTIC GAS 





The equation of state for a completely degenerate non-relativistic 
gas (i.e., one in which pp < mec’) is fairly straightforward. When 
degenerate, the number of electrons with momenta between p and 
p+dp is 
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To calculate the pressure (the flux of momenta through a unit 
surface), you need to know two things: 


1) The flux of electrons hitting the surface, and 
2) How much momentum is transferred by each electron. 


The former quantity can be calculated in a manner similar to dif- 
fusion. Recall that, for diffusion, we said that the flux of particles 
with velocity v diffusing across a boundary from one direction is 
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where l is the particle mean free path and v is the particle veloc- 
ity. If the mean free path is small, the total number of electrons 
crossing the boundary from both directions is 
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If each electron carries an amount of momentum, p, then the 
pressure at the boundary from electrons with velocity v is 
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and the total pressure is 
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For a non-relativistic electron gas, v = p/Me, so for complete 
degeneracy 
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If we then substitute for py using (7.3.3) 
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Note that there is no temperature dependence in this equation. 
The distribution of electrons in phase space is entirely determined 
by the Fermi exclusion principle, leaving pressure proportional 
only to density. 








Since there must be a smooth transition between the degenerate 
and non-degenerate states, we can estimate the density where this 
occurs by simply equating the two expressions from pressure 


h2 2/3 5/3 
Frais = Lopp Soe ak ( à ) (2) 
HUHMa 20MeMa NMa He 


Degeneracy therefore becomes important where 
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COMPLETELY DEGENERATE, RELATIVISTIC GAS 


If the density is extremely high, the velocities that the electrons 
must have to obey the Pauli exclusion rule will be relativistic. In 
this case, we must substitute for velocity using 
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The integral for pressure is then slightly more complicated 


3 h’ 1 + p?/(m2c? 


(7.3.9) 





If we let € = p/(mec), and x = py/(mec), then the integral can 
be re-written as 
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Following Chandrasekhar’s Principles of Stellar Structure, if we 
let 
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Note that if the gas is extremely relativisitic, py >> mec, thus 


x >> 1, and f(x) — 22+. In this limiting case 
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Substituting density for pr using (7.3.6), we get 
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Again, the pressure is independent of temperture, but this time, 


the exponent is 4/3. 


The transition point between the relativisitic and non-relativistic 
cases occurs when pf ~ mec. If we substitute mec into (7.3.6), we 


find that this occurs at 
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pP > ps ee = 9.7 x 10° g cm~’ 
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PARTIAL DEGENERACY 


In many cases, when both the temperature and density are high, 
partial degeneracy occurs. When this happens, the electron dis- 
tribution in phase space is neither Maxwellian 
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nor degenerate 
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where w is the degeneracy parameter, 
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and EF is the kinetic energy of the electron 
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Note that this equation does have the correct form. If ne is small, 
then Y <0, and ff (p) — f~ (p). On the other hand, if ne is 
large, e@/*T-¥ < 0, and f? (p) — fP (p). 


Note also that 7 depends only on n./T?/?. 





With this general expression for the electron phase-space occupa- 
tion, the expressions for ne, Pe, and ue are 
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where in the non-relativisitic case, 
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For the non-relativistic case, if we let n = p?/2m,kT, then these 
integrals become 
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All the above integrals have the form 
F (p) = J i S dn (7.3.25) 


These are called Fermi-Dirac integrals, and they are tabulated in 
McDougall & Stoner 1939, Phil. Trans. R. Soc. London, 237, 67 
(so you don’t need to numerically integrate them each time). 


The General Equation of State 





We get the general equation of state for stellar matter by summing 
the expressions for radiation pressure, ion pressure, and electron 
pressure. ‘Thus 
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Thus, we have implicit formulae for the equation of state. Given 
p and T, you can solve for via (7.3.16), and then P. Similarly, 
the equation for each particle’s internal energy is 
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Unfortunately, there is no way to use these equations to generate 
analytical expressions for the thermodynamic quantities cp, cy, 
a, ô, and Vaa. Unless one of the terms dominates, these quantities 
must be computed numerically. 


